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ABSTRACT 



Multi-zone chemical evolution models (CEMs), differing in the nucleosynthesis prescriptions (yields) and prescriptions of star for- 
mation, have been computed for the Milky Way. All models fit the observed 0/H and Fe/H gradients well and reproduce the main 
characteristics of the gas distribution, but they are also designed to do so. For the C/H gradient the results are inconclusive with 
regards to yields and star formation. The C/Fe and 0/Fe vs. Fe/H, as well as C/0 vs. 0/H trends predicted by the models for the solar 
neighbourhood zone were compared with stellar abundances from the literature. For 0/Fe vs. Fe/H all models fit the data, but for C/O 
vs. O/H, only models with increased carbon yields for zero-metallicity stars or an evolving initial mass function provide good fits. 
Furthermore, a steep star formation threshold in the disc can be ruled out since it predicts a steep fall-off in all abundance gradients 
beyond a certain galactocentric distance (~ 13 kpc) and cannot explain the possible flattening of the C/H and Fe/H gradients in the 
outer disc seen in observations. Since in the best-fit models the enrichment scenario is such that carbon is primarily produced in 
low-mass stars, it is suggested that in every environment where the peak of star formation happened a few Gyr back in time, winds of 
carbon-stars are responsible for most of the carbon enrichment. However, a significant contribution by zero-metallicity stars, espe- 
cially at very early stages, and by winds of high-mass stars, which are increasing in strength with metallicity, cannot be ruled out by 
the CEMs presented here. In the solar neighbourhood, as much as 80%, or as little as 40% of the carbon may have been injected to the 
interstellar medium by low- and intermediate-mass stars. The stellar origin of carbon remains an open question, although production 
in low- and intermediate-mass stars appears to be the simplest explanation of observed carbon abundance trends. 
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1. Introduction 

Carbon is one of the most common elements, but we know sur- 
prisingly little about its origin. What we do know, however, is 
that it is ubiquitous throughout the Universe, and can be found 
in just about any astrophysical environment. We also know that 
life, as we know it, requires the existence of carbon, nitrogen, 
oxygen and a few other elements. Understanding the origin of 
carbon may therefore tell us something about the probability of 
finding carbon-based life elsewhere in the Galaxy, i.e., beyond 
the solar neighbourhood. 

The stellar origin of carbo n is mainly due to the Triple- 
Alpha reaction ('Salpet er 1952h but this reaction may occur in 
various types of stars. Carbon Stars (C-stars) have been recog- 
nised as a class of astronomical object for more than a cen- 
tury and have several times been suggested as the main carbon 
sources in the Universe. Already in the work by Burbidge et al. 
(1957) it was suggested that carbon was provided by mass-loss 
from red giants and supergiants. Later Dearborn (1978) sug- 
gested that low-mass stars may be a significant source of car- 
bon based on abundance determinations in planetary nebulae. 
Recent theoretical work on stellar evolution of low and interme- 
diate mass (LIM) stars also points to low-mass stars as being 
significant producers of carbon, although the quantative results 



signincant proaucers or carbon, aitnougn tiie quantative results 
may differ (fa arigo 200j yizzard et al. 2004tlGavilan et al. 2001 
iKarakas & Lattanzio 200% ! 
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Models of chemical evolution (CEMs) are in general in good 
agreement with observed abundances if a delayed carbon re- 
lease from LIM-stars is assumed. Timmes et al. (1995) used 
the nucleosynthetic yields by Woosley & Weaver (1995, from 
hereon cited as WW95) and Renzini & Voli (1981), which 
led to a very significant contribution of carbon from LIM- 
stars to the Galactic Disc, and more recent work (using other 
sets of yields for LIM-star s) have led to q t iite similar results 
(IChiappini et al. 2C)03: .Akerman et al. 2004t ICarigi et al. 20051: 
iGavilan et al. 2005h . However. Maeder (1992) argued that ra- 
diatively driven winds from high-mass (HM) stars should pro- 
vide huge amounts of helium and carbon. In such case, these 
stars would be the main contributors, and in CEMs the role 
of LIM-stars would have to be much less significant to avoid 
over-production of carbon compared with observed abundances. 
Garnett et al. (1995) observed that the C/O-ratio increased with 
increasing O/H in dwarf irregular galaxies, which they inter- 
preted as being consistent with carbon being produced in HM- 
stars with metallicity-dependent yields, as in the models by 
Maeder (1992) and Portinari et al. (1998). Following that idea, 
Gustafsson et al. (1999) argued that the rising C/O-trend with 
metallicity that they found in Galactic-disc stars was the result 
of carbon being produced in HM-stars rather than LIM-stars. 

Recent observations have revealed a declining trend with 
increasing metallicity for the C/O ratio in the solar neighbour- 
hood at early times (Akerman et al. 2004; Fabbian et al. 2009aj). 
C/O vs. O/H shows a negative slope roughly until the an- 
ticipated onset of disc formation, i.e., during the first billion 
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Hubble pi'ofile 




= 15.0 kpc 



Fig. 1. The gassian surface density profile adopted for the 
halo/thick disc in this paper compared to a Hubble (1923) pro- 
file, normalised to the density at the Sun's galactocentric dis- 
tance (marked by a vertical dashed Une in the figure). 

years of Galactic evo l ution, which is in d isagreement (see, e.g. 
IChiappini et al. 2003t lOavilan et al. 2005h with the predictions 
of CEMs that do not include any modifications of the stan- 
dard carbon and oxygen yields. Such C/O discrepancy is ob- 
viously connected to the first generations of stars in the Milky 
Way, and may therefore have common origin with our still in- 
complete understanding of early chemical evolution. For in- 
stance, the underabundance of essentially metal-free LIM-stars 
in the halo, which is often claimed to be a result of a top- 
heavy initial mass function (IMF) at early times (s ee, e.g., 
lAbel et al. 200l iTumUnson 2005 iKarlsson et al. 2008L and ref- 
erences therein), is one such problem. If the IMF has evolved 
from being initially top-heavy, to the form that is observed in 
the solar neighbourhood today, it is possible that this may affect 
the evolution of the C/O ratio as well. Chiappini et al. (2000) 
have considered an evolving IMF, but concluded that it did not 
improve the agreement with observational constraints. However, 
their study did not focus on the early chemical evolution and 
abundance trends at very low metallicity. 

In this paper, the origin of carbon is investigated once again, 
using a set of multi-zone CEMs for the Milky Way. Several nu- 
cleosynthetic prescriptions are considered, as well as the effects 
of an evolving IMF. 

2. Chemical evolution model 

2.1. Galaxy Formation 

It is usually assumed that the Galaxy was formed through 
baryonic infall, or more precisely, by accretion of pristine gas 
(hydrogen and heliu m), and that the rate of accretion follow s 
an exponential decay (iLacev & Fall 19851 iTimmes et al. 19951) . 
Furthermore, it is rather well estabUshed that the halo/thick disc 
and the thin disc components of the Galaxy were assembled 
on different time-scales and, perhaps, with some separation in 



time, as suggested by Chiappini et al. (1997). The latter scenario, 
known as the two-infall model, consists of two infall episodes 
where the disc formation starts after some time rmax (which is 
when the rate of accretion reaches its maximum), i.e., the total 
rate of accretion is 



£inf.(nO 



2h(r, f) if f < Tmax 

ih(n + 2d(r, f) if f > r„„ 



where, for the halo/thick disc, 

2h(r,f) = 
and for the (thin) disc, 

^d('-,ro) [, [ (^0 - Tmax) 



2h(r,fo) 

-i 1 - exp 
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2;d(r,f) = 



^d(r) 



1 - exp 



X exp 



(f- 



(3) 



where r is the infall time scale and t - fo is the age of the 
Galaxy. Here, the value used for fo is 13.5 Gyr, which can be 
regarded as an upper limit. To simulate the inside-out forma- 
tion of a galactic disc, the infall time scale is assumed to have 
a linear radial d e pendence (for further details, see Chiosi 
iMatteucci & Francois 19891 and references therein), 



Tci(r) = max 



O,r0-T„l 



Re 



(4) 



where is the infall time scale in the solar neighbourhood, tq 
is an arbitrary constant and Rq - 8.5 kpc is the galactocentric 
distance of the Sun. For the halo/thick-disc phase it is assumed 
that Th = 1.0 Gyr at all galactocentric distances. 

The present-day total baryon density in the solar neighbour- 
hood is assumed to be 6 0Mo pc^ (which is consi stent with the 
results obtained by, e.g.. lHolmberg & Flvnn 2004h and the final 
baryonic (thin) disc is assumed to follow an exponential distri- 
bution, 

Sd(r,fo) = So(fo)expj-^j, (5) 

where So is the present-day central surface density and is the 
disc scale length. The value of the latter parameter is not very 
well constrained for the Milky Way. In this study it is assumed 
that is rather short, which is supported by stellar statistics 
(Ruphy et al. 1996; Porcel et al. 1998). The final density Eh of 
baryonic matter of the halo/thick-disc component in the solar 
neighbourhood is of the order of lOM© pc"^. The present-day 
density of halo stars in the s olar neighbourho od has been esti- 
mated to 5.7 ■ 10 "^Mq kpc"^ (iPreston et al. 1991 ). which coitc- 
sponds to a surface density of 1 - 2Mq pc~^ and for the thick 
disc stars the local surf ace density is about 3.5% of that of 
the thin disc (tOhia 200 lb . But remains difficult to constrain 
since the corresponding gas fraction is currently unknown. The 
halo/thick-disc surface density is modelled by a gaussian. 



Y.h(r, f) = Sh(^0, fo) exp 



R, 



•1 \ 



(6) 



where R^ is defined as above and rh is the scale length. The inner 
part of this profile is si milar to a Hubble ( 1923) profile, used in 
some other studies (e.g., Renda et al. 2005), while the outer part 
is gradually steeper (see Fig. [Hi, as suggested by very deep star 
counts in the halo (see Hel mi 20081 and references therein). For 
further details, see Table [1] 
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Table 1. Parameters for the CEMs. 



Model 


2h(^o,ft)) 


''h 


"^diJ^Q, to) 


rd 




Vd 




£ 








Mo pc"- 


kpc 


Ma pc"^ 


kpc 


Gyr' 


Gyr-' 


Gyr-' 




Me 


Me 


Al 


10.0 


15.0 


50.0 


2.1 


0.170 


0.170 


0.0425 


1.5 


112.0 


0.32 


Bl 


10.0 


15.0 


50.0 


2.1 


0.138 


0.138 


0.0345 


1.5 


100.0 


0.40 


CI 


10.0 


15.0 


50.0 


2.1 


0.163 


0.163 


0.0408 


1.5 


113.0 


0.32 


Dl 


10.0 


15.0 


50.0 


2.1 


0.166 


0.166 


0.0415 


1.5 


109.0 


0.32 


El 


10.0 


15.0 


50.0 


2.1 


0.176 


0.176 


0.0440 


1.5 






A2 


15.0 


11.5 


45.0 


2.1 


0.535 


0.535 


0.0000 


1.0 


114.0 


0.30 


B2 


15.0 


11.5 


45.0 


2.1 


0.404 


0.404 


0.0000 


1.0 


100.0 


0.40 


C2 


15.0 


11.5 


45.0 


2.1 


0.535 


0.535 


0.0000 


1.0 


114.0 


0.30 


D2 


15.0 


11.5 


45.0 


2.1 


0.535 


0.535 


0.0000 


1.0 


110.0 


0.30 



2.2. Star formation 



2.3. The stellar initial mass function 



For the halo/thick disc, the star-formation rate is prescribed by a 
modified Schmidt-law of the form 



i*(r,f) = VhS(r,f) 



'^^asir, t) 



S(r,f) 



(7) 



where Vh is the star-formation efficiency expressed in Gyr-' . For 
the halo/thick disc phase it is assumed that Vh is constant, while 
for the disc phase the star-formation efficiency is assumed to be 
related to the angul ar frequency of the disc dWvse & Silk 1989i: 
iBoissier et al. 2003h . i.e.. 



S*(r,f) = v(r,f)2(/?Q,fo) 
and 



to) 



(8) 



v(r, f) = 



Qo(r) 



if Saas > 



if Xgas 



(9) 



where Sc;as is the gas density, is the critical density, Vd, are 
constants and Qq is the present-day (f = fo) angular frequency 
of the disc (rotation data were taken from Sofue et al. 1999). To 
avoid unphysical "kinks" when switching between the two star 
formation regimes, a "smooth step function", defined as. 



0(r,t) = 



1 



1 H- exp[-2^(2gas -Sc)] 



(10) 



is used instead of an actual conditional expression. This func- 
tion will approach a regular Heaviside step function as ^ —> oo, 
but here it is assumed that ^ = 1. The critical density Ec may 
change along the disc, but for sim plicity it is assu med that 

= 7.OM0 pc-2 over the whole disc (iKennicutt 19891) . Two dif- 
ferent star formation prescriptions are considered in this paper: 
(1) a Schmidt law as above with s - 0.5 (cf. Fuchs et al. 2009) 
above the critical gas density 2c, and below Sc star formation 
proceeds with lower efficiency according to a linear Schmidt- 
law ii, = VjSgas, where - 0.25vd (models Al to El in Table 
[TJ, and (2) a Schmidt law with e = and a strict threshold, i.e., 

= ^ 1* = below the critical density (models A2 to D2 in 
Table [1). The first case will be referred to as star formation law 
of type 1, and the second case as type 2. 



It is assumed that stars are formed according to a stellar IMF 
that may, or may not, be time-dependent. In its simplest form, 
the IMF is just a power-law with sharp cut-offs at the lower and 
upper ends. It is well established, however, that the IMF turns 
over at low masses and is probably truncated at the high-mass 
end (see Fig. |2] and references therein). Hence, an IMF of the 
form 



(ot, t) — (poit) m 



-(\+x) 



exp 



OTc(0 



'«u(f) 



(11) 



is adopted, where m^, my are the masses defining the low-mass 
turn over and the high-mass truncation of the IMF, respectively. 
The two parameters m^, may be regarded as functions of time, 
thus allowing for an evolving IMF. The constant (p^, which is 
obtained by the normalisation condition 



Jo 



m cf)(m, t) dm — 1 , 



is then also a function of time. More precisely, 
1 /^'(f) 



<^o(0 



x>\. 



2 K,^i{^l) 

where is a dimensionless variable defined as 



^l{t) = 



nicit) 



\ tnu(t) ' 



(12) 



(13) 



(14) 



and K„ is the modified Bessel function of the second kind and 
order n. The mean stellar mass of this IMF is 



<m)(f) = niMnit) 



x>l. 



KM ' 

while the most probable mass is given by 
1 



tna{t) - -xOTu(0 



■^l+4x^fi(t)^ 



1 



(15) 



(16) 



Both the mean and most probable masses are also functions of 
time in the general case. However, in the special case when 
nicit) oc m^it), or - const., the situation is somewhat simpler 
First of all, 0o becomes constant over time and the IMF does 
not need to be renormalised at every time step. This simplifying 
assumption will be used throughout the following. 

The time-evolution of the parameters and cannot be 
completely arbitrary. They are quite likely related to some char- 
acteristic mass-scale related to the physical origin of the IMF. 
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Fig. 2. Comparison between the IMF 0(ot) used in this paper and 
some IMFs frequently used in the Hterature. The stellar mass m 
is given in solar masses. 



Larson (1995; 1996; 1998) pointed out that there should be a 
connection between the turn-over mass (or the characteristic 
stellar mass) and some fundamental mass-scal e in the star for- 
mation process, such as the Jeans mass (Jea ns 1902h . Recent 
observational evidence for an evolving IMF suggests charac- 
teristic stellar mass es which are consistent with such a picture 
dvan Dokkum 20081) . The thermal Jeans mass is given by the 
pressure and temperature of the collapsing cloud, i.e., mj oc 
j2p-\l2^ In a self-gravitating cloud, the pressure-gravity balance 
is such that the Jeans mass can also be expressed in terms of the 
gas surface density 2 as mj cc T'^Yt'^ (see, e.g., Larson 1985). If 
the ISM is isothermal, the Jeans mass is inversely proportional to 
the local surface density of the gas and if cooling is inefficient, it 
may be close to constant in time. In the following we assume that 
the gas temperature J is a function of the gas density S alone, 
i.e., a polytropic equation of state. Assuming then that and 
are related to the characteristic Jeans mass of the ISM, it is rea- 
sonable to parametrise these quantities in terms of the local gas 
mass density of the Galactic disc (.Elmegreen 1999) . In particu- 
lar, a power-law form, m^{t) oc ;«„(?) oc 2gfs(f), with p e [0, 1], 
gives an adequate description of how and change during 
the evolution of the Galaxy (see Fig. [3] showing how the IMF 
changes with the gas density). 

In the present study, two types of IMFs according to Eq. ( fTTT i 
are considered: 

(a) A non-evolving IMF of the form 



^{m) - 00 ni exp 



- 


'OTc 


m \ 




^ m 





(17) 



where the upper mass-cut is 100 - I2OM0 in order to avoid 
over-production of oxygen, and the turn-over mass is 
0.30 - O.4OM0. The power-law i ndex x w as chosen to be 
steeper than the canonical (Salp eter 1955h value, i.e., x - 
1.80 instead of x = 1.35, in order to reproduce the proper- 



Fig. 3. Evolution of the IMF with the gas density in model El. 
Meaning of the axes as in Fig.|2l 



ties and abundances of the solar neighbourhood (see, e.g., 
Chiappini et al. 1997). 
(b) A time-dependent IMF of the form given in Eq. (fTTT l. where 



'Wc(0 



and 



OTu(f) 



»Jc(fo) 



. »Jc(fo) 



mu(fo) 



, OTu(fo) 



-gas 



4M0PC- 



4Mopc- 



-0.6 



if Sgas < 4M0PC" 



if Sgas > 4M0PC" 



(18) 



-0.6 



if Sgas < 4Mopc 



if Igas > 4Mopc" 



(19) 



In the equations above, mc(fo) and mu(fo) denote the present- 
day values for the local-disc IMF. Note that when Sgas > 
AMq pc"^, this IMF is nearly identical to that of case (a). 

A comparison between the IMF adopted here and a few IMFs 
commonly used in CEMs is shown in Fig.|2l Note that case (b) 
above will lead to a larger fraction of very low-mass stars and 
substellar objects, compared to, e.g., the IMF by Scalo (1986), 
which means that a larger fraction of metals will be locked-up in 
stellar objects that do not release any significant amounts of gas 
back to the ISM within a Hubble time. 

In order to include the contribution from supernovae type la 
(SNIa) events, which is commonly assumed to be the result of 
gas accretion onto a white dwarf from its companion (presum- 
ably a red giant) in a binary system, the formalism introduced 
by Greggio & Renzini (1983) and Matteucci & Greggio (1986) 
is used. In this formalism, the equation describing the evolution 
of an element / can be written as 

ot Z(r, f) 



f 

J16 



^i{r, t - Tm) i/>(m) i-^{r, t - t^) dm + 
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Fig. 4. Predicted metallicity distribution for the solar neighbour- 
hood G-dwarfs according to models Al-El convolved with a 
gaussian with cr = 0.15 to simulate observational errors. 
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Fig. 5. Cumulative fraction of carbon due to LIM-stars as a func- 
tion of time at the solar vicinity, according to model Al-El in 
this paper 



r Xi{r,t-T,„)(/>(m)i^(r,t-T,„)dm+ (20) 
Jo 

X16 
Xi{r, t - T,n) (/>im) t*(r, t - t„) dm + 

plk ^1/2 

T] I I (f(ju)Xi(r,t -Tm)(p(m)ti,{r,t -T„)dfidm, 

J 3 Jfimin 

where Tm is the lifetime of a star of mass m (as obtained from the 
stellar evolutionary tracks by Schaller et al. 1992 and WW95), 
(f>{m) is the IMF described above, X, is the production matrix for 
an element / as defined by Talbot & Arnett (1971, 1973) and rj 
is the fraction of (binary) low and intermediate mass stars un- 
dergoing a SNIa event. In the equation above, it is implicitly 
assumed that SNIa events originate from binary systems more 
massive than 3Mq and with a maximum mass of 16Mq. The bi- 
nary distribution function (fi/u) is a function of the ratio of the 
secondary to the combined mass of the binary system fi, which 
is normalised on the interval yU E [0, 0.5]. 

3. Abundance data 

For the Galactic HII regions, the oxygen and carbon abundances 
and the distances relative to the Sun, were taken from Esteban 
et al. (2002; 2005). For the extragalactic HII regions, data were 
taken from Bergvall (1985), Garnett et al. (1995; 1997; 1999), 
Kobulnicky et al. (1997) and Izotov & Thuan (1999). In all cases 
the spectra are corrected for extinction and underlying absorp- 
tion in the Balmer lines. The carbon abundances by Esteban et 
al. (2005) are also corrected for the presence of a carbon dust 
component. 

Stellar abundance data were compiled from several 
different studies. For the solar neighbourhood, the data 
were compiled from several authors (Bensby & Feltzing 2006; 
lEcuvillonet al. 2004t iFabbian et al~2009a; .Gratton et al. 2000c 



[Gustafss on et al. 1 9991; ll sraelian et al. 2 004 iJonsell et al. 20051; 
Spite et al. 2005), while the Cepheid abundances in Fig. |7] were 
taken from the works by Andrievski et al. (2002a; 2002b; 2002c; 
2004) and Luck et al. (2003; 2006). The observed G-dwarf dis- 
tributions used in this paper were taken from Edvardsson et al. 
(1993) and Wyse & Gilmore (1995). 

Results for HI and H2 (CO) were taken from Dame (1993). 
To obtain a consistent and, possibly, more correct set of data for 
the hydrogen distributions, the conversion from CO flux to H2 
surface density has been redone using the oxygen-dependent cal- 
ibration by Wilson (1995), assuming the 0/H-gradient derived 
by Esteban et al. (2002; 2005). 

4. Results and discussion 

Using the framework described in Sect. |2] in association with a 
variety of nucleosynthetic prescriptions, we computed the fol- 
lowing five types of CEMs for the Milky Way: 

(A) HM-star yields by WW95 and LIM-star yields by Gavilan 
et al. (2005) and Gavilan et al. (2007), the latter being for 
low-metallicity. 

(B) HM-star yields by Portinari et al. (1998) and LIM-star yields 
from van den Hoek & Groenewegen (1997). 

(C) As in model A, except that for HM-stars the carbon yields at 
zero metallicity (Z = 0) are increased by 50%. 

(D) As in model C, except that the LIM-star carbon yields below 
Z - 0.01 are increased by 50% in an attempt to reproduce the 
highest carbon abundances derived for objects having metal- 
licities around solar. 

(E) Yields as in model A but now with an evolving IMF accord- 
ing to Eq. ( fTSl l. 

Case A represents a case where the LIM-stars are producing 
most of the carbon, while case B represents a case where the 
carbon is to a large extent produced in HM-stars, similar to the 
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O Gustafsson et al. (1998) 

+ Ecuvillon et al. (2004) 

>K Fabbian et al. (2009) 

OBensby & Feltzing (2006) 

AJonsell et al. (2005) 

□ Spite (2005) 

X Israelian et al. (2004) 

• Gratton et al. (2000) 

end type of star formation law (A2-D2) are not 
siiown, but yield very similar results. The dis- 
continuity in the model tracks seen in the car- 
bon abundance plot around [Fe/H] = -1 is due 
to the onset of disc formation, i.e., the second 
infall episode. 



Model Al 

Model Bl Fig. 6. Stellar abundances of carbon, nitrogen 

Model CI and oxygen, relative to iron in the solar neigh- 

Model D 1 bourhood. All lines (models) and symbols (ob- 

Model El servations) are explained in the figure. Data 

from Fabbian et al. (2009a) is that derived ac- 
counting for H-collisions. The thin dashed lines 
marks the solar values. Models with the sec- 
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Table 2. Predicted and observed quantities of the solar neighbourhood. 2* and lint denote star-formation rate and infall rate, 
respectively. 



Quantity 


Unit 


Al 


Bl 


CI 


Dl 


El 


A2 


B2 


C2 


D2 


Observed 


SN I rate 


pc-^ Gyr' 


1.11 


1.44 


1.15 


1.12 


1.09 


1.00 


1.54 


1.00 


1.00 




SN II rate 


10-3 pc-2 Gyr-' 


3.81 


3.75 


3.87 


3.89 


3.69 


3.17 


3.82 


3.17 


3.20 




SN II/SN I 




3.43 


2.60 


3.37 


3.47 


3.39 


3.17 


2.48 


3.17 


3.20 


~ 3.1 




Mo pc-- Gyr-' 


3.25 


2.82 


3.30 


3.33 


3.17 


2.82 


2.89 


2.82 


2.87 


2-5 


£gas(^0> to) 


Mo pc-2 


7.81 


7.98 


7.95 


7.93 


7.68 


7.44 


8.04 


7.44 


7.47 


8±5 


^mf(Ro, to) 


Mo pc-2 Gyr-' 


1.44 


1.44 


1.44 


1.44 


1.44 


1.30 


1.30 


1.30 


1.30 


0.5-5 



Table 3. Solar abundances by mass predicted for 4.56 years ago by the models and compared with the current solar composition 
derived from observations. 





Al 


Bl 


CI 


Dl 


El 




A2 


B2 




C2 


D2 




Asplund et al. (2005) 


H 


0.736 


0.729 


0.736 


0.736 


0.736 




0.736 


0.729 




0.736 


0.736 




0.739 


He 


0.253 


0.260 


0.253 


0.253 


0.253 




0.253 


0.260 




0.253 


0.253 




0.249 


C 


2.18 10-3 


2.77 10-3 


2.18 10-3 


2.85 10-3 


2.18 10" 


-3 


2.18 10-3 


2.69 10 


-3 


2.20 10-3 


2.86 10- 


-3 


2.18 10-3 


N 


9.32 10-* 


7.54 10-" 


9.24 10-" 


9.43 10-" 


9.32 10- 


-4 


9.22 10-" 


7.38 10- 


-4 


9.22 10-" 


9.47 10- 


-4 


6.23 10-" 





5.41 10-3 


5.41 10-3 


5.41 10-3 


5.41 10-3 


5.41 10 


-3 


5.41 10-3 


5.41 10- 


-3 


5.41 10-3 


5.41 10- 


-3 


5.41 10-3 


Fe 


1.17 10-3 


1.17 10-3 


1.17 10-3 


1.1710-3 


1.1710 


-3 


1.17 10-3 


1.17 10- 


-3 


1.17 10-3 


1.17 10- 


-3 


1.17 10-3 



scenario suggested by Gustafsson et al. (1999). Case C is an at- 
tempt to explain the C/O trend obtained by Akerman et al. (2004) 
and Fabbian et al. (2009a) from observational data of stars at low 
metallicity. The modified LIM-star yields in Case D aims at ex- 
plaining the super-s olar carbon abundance s observed in the so- 
lar neighbourhoo d ("Gustafsson et al. 1999'; 'Esteban etaLjOOl 
ICarigi et al. 2005t iBensbv & F eltzing 2006). Finally, case E is 
introduced as an alternative solution to the early C/O evolution, 
where instead of modifying the yields for the first generation of 
stars, a top-heavy IMF is used at low gas densities. In all mod- 
els the iron yields of WW95 at Z = have been reduced by a 
factor of five in order to obtain O/Fe-ratios consistent with those 
found in very metal-poor stars (see also the discussion in Sect. 
14.4b . That the yields of WW95 may lead to an over-production 
of iron, at low metallicity in particular, was pointed out already 
by Timmes et al. (1995). 

The model parameters were adjusted to achieve good agree- 
ment with the abundance data described in Sect. 3, i.e., either at- 
tempting to reproduce the radial trends of elemental abundances 
(star formation type 1, referred to as models Al-E l, see Fig.lTT l 
or of the molecular hydrogen density as derived by (iDame 1993h 
from observations of the outer parts of the Galaxy (star forma- 
tion type 2, referred to as models A2-D2, see Fig.[8]l. Other prop- 
erties, such as the solar abundance pattern, the observed G-dwarf 
distribution and basic solar neighbourhood quantities were also 
used as constraints in both cases (see Table |2l Tableland Fig. 
II. 

4.1. Solar neighbourhood trends 

For the chemical evolution of the solar neighbourhood (defined, 
in the models presented here, as the evolution at galactocentric 
radius R - 8.5 kpc) prescription of type 1 and 2 for star forma- 
tion provide very similar results. Hence, only models Al-El will 
be discussed here. Model Al provides satisfying results for the 
O/Fe trend and also a quite good agreement with the observed 
C/Fe trend, except at disc-like metallicities ([Fe/H]> -1), where 
observations suggest a higher C/Fe (see Fig.|6]l. However, model 
Al fails to explain the C/O trend shown by metal-poor stars (Fig. 



|9]l. Model Bl, in which a considerably larger fraction (compared 
to model Al, see Fig.|5]l of carbon comes from HM-stars seems 
to slightly over-produce carbon relative to iron, compared to cor- 
responding ratios derived from observations for the solar neigh- 
bourhood and, just as model Al, it cannot reproduce the C/O 
trend at low metallicity. 

In models CI and Dl, the carbon production at zero metal- 
licity is significantly increased, which may provide a solution 
to the carbon-enrichment problem at low metallicities, as sug- 
gested by Akerman et al. (2004) and Carigi et al. (2005). The 
Z = yields for HM-stars by Chieffi & Limongi (2004), provide 
a C/O-ratio that is consistent with that derived from observa- 
tions at very low metallicity, but it is roughly the same also at 
moderately low metallicities where observations instead suggest 
that C/O becomes much lower. If the carbon yields were indeed 
higher than according to WW95, but only at Z = 0, it would 
pose an interesting problem for nucleosynthesis modelling: why 
do massive, essentially metal-free, stars produce much more car- 
bon than more metal-rich ones of similarly high mass? 

In the present paper, effects of inhomogeneities and con- 
tributions from pair-instability supernovae during the earliest 
evolution of the Galaxy are not studied, although the lat- 
ter may still play an importan t role also for the C /O-ratio 
(cf. Karl sson & Gustafs son 2005: Karlsson et al. 200 8*). The to- 
tal carbon yield may thus be significantly higher at early times, 
as in models CI and Dl. Despite the high carbon yields at Z = 0, 
also models CI and Dl are reasonably consistent with the ob- 
served C/Fe vs. Fe/H-trend at low metallicity (again, see Fig.|6]l. 
Model Dl, however, is the only model that seems to explain the 
highest C/Fe-ratios found at the high metallicity end. 

The solar abundance pattern (taken to be the composition 
of the ISM at R - Rq, 4.56 Gyr ago) is well reproduced by 
all models - except for an over-production of nitrogen in gen- 
eral and for models Dl and D2, where the carbon abundance is 
about 30% too high (see Table O - but it should be noted that 
all model parameters are calibrated such that solar abundances 
and the peak of the G-dwarf distribution (Fig.|4ll are reproduced. 
For disc-like metallicities, the apparent over-production of nitro- 
gen by the models seen in Fig.|6]may be a result of the LIM-star 
yields chosen for this study in combination with the prescription 
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Fig. 7. Predicted and observed radial trends for O/H, C/H (upper panels), Fe/H and the surface density of hydrogen (lower panels) 
in the Milky Way for models computed using the star formation prescription of type 1 . Solid circles show abundances in HII regions, 
while empty circles show abundances in Cepheids. The vertical and horizontal thin, dashed lines indicate the solar values according 
to Asplund et al. (2005). 



for star formation. Both van den Hoek & Groenewegen (1997) 
and Gavilan et al. (2005) predict high nitrogen yields. In the low- 
metallicity end, however, the models seem to underproduce ni- 
trogen, although the high nitrogen abundances derived from ob- 
servations of some hal o stars are probably the exception rather 
than the rule ( Israelian et al. 2004 . 

Model El, in which an evolving IMF was used, provides 
an alternative solution to the early carbon-enrichment problem. 
Using the yields by WW95, without modifications, it reproduces 
the C/O vs. O/H trend and simultaneously the C/Fe vs. Fe/H 
at low metallicity fairly well. The reason why an evolving IMF 
can solve C/O vs. O/H problem using the yields by WW95, is 
that the C/O-ratio of the ejecta is much higher for stellar masses 
around = 10 - 15M0 than for more massive stars. Hence, if 
the IMF is peaking at intermediate stellar masses rather than at 
subsolar masses, the IMF weighted C/O-ratio is roughly 0.5 dex 
higher than it would be otherwise. Whether an evolving IMF is 
more likely than a scenario with high carbon yields at Z = 
is matter of debate, but several recent studies suggest that in 
the early Universe, star formation took place according to an 
IMF that was more top-heavy than at present time (see, e.g., 
iTumlinson 20061 : [Pave 20081 : Ivan Dokkum 20081 and references 
therein). 

4.2. Abundance gradients 

The Fe/H-gradients predicted by the models all agree quite well 
with the observed Fe/H-gradient obtained from Cepheids (see 
Fig. IT). A similar radial trend was also found by Nordstrom et 



al. (2004) for young stars at galactocentric distances between 6 
and 10 kpc. 

The results regarding the C/H-gradien t are inconclusive . 
The ISM abundances (Est eban et al. 2005t ICarigi et al. 20051) 
are quite high and in reasonable agreement with all models, al- 
though models CI and C2 fit these data better (see Fig. |7] and 
IS). On the other hand, the Cepheid carbon abundances are sig- 
nificantly lower, which better matches, at solar galactocentric 
distance, the current solar abundance of carbon. However, the 
ranges of abundances seen in in the solar neighbourhood are 
rather wide and the Sun may therefore not represent the typi- 
cal composition of the ISM in the solar neighbourhood at the 
time of formation of the Sun. This suggests that calibrating the 
parameters of the CEMs against the current solar abundance pat- 
tern may cause the present-day C/H predicted by these mod- 
els to end up above (or below) that measured in young stars, 
which one may assume reflects the C/H of the ISM. It might 
be that the solar abundance pattern should not at all be used as 
a constraint on CEMs for the epoch when the Sun was formed, 
since its chemical composition (relative to iron) may depart from 
those of most solar-t ype stars of s imilar ages and orbits (for fur- 
ther discussion, see iGustafsson 2008. an d references therein) , 
in fact even from those of solar twins ( Melendez et al. 20091 : 
iMelendez et al. 2009; Nordlund 2009). In any case, the abun- 
dance difference between Cepheids and Hll-regions still lacks 
a proper explanation. 

Models A2-D2 were designe d to fit the h ydrogen distribu- 
tion (according to the results of iDame 1993h as well as pos- 
sible. The good agreement is mainly a result of the assump- 
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Fig. 8. Same as Fig.|7] but for models computed using the star formation prescription of type 2. 



tion that there exists a true star-formation threshold in the 
disc, below which star formation will essentially cease com- 
pletely (in princip le identic al to the assum ption made by, e.g., 
IChiappini et al. 19 97; Chia ppini et al. 2001h . However, the fit to 
the hydrogen distribution is obtained at the expense of a good 
fit to the abundance gradients derived from observations in the 
outer parts of the Galaxy. As shown in Fig. [8] the threshold 
causes the abundance gradients predicted by the models to fall- 
off" steeply at galactocentric distances beyond 13 kpc, which is 
the part of the disc where the gas density in the model never 
reaches the threshold value. On the basis of this disagreement 
with C, O and Fe abundance gradients derived from observa- 
tions, it is therefore not likely that star formation will actually 
cease completely when the gas density drops below the critical 
density. 

The somewhat poorer agreement between models A 1 -E 1 and 
the observed hydrogen distribution should not be seen as a ma- 
jor problem. In the present study, the possibility of having a slow 
radial gas flow has not been studied. Such a flow has not been 
observed, but even a very slow radial gas flow can alter the hy- 
drogen distribution quite significantly, provided enough time is 
available for the process. Hence, the hypothesis presented here 
via models Al-El (i.e., that there are two modes of star forma- 
tion rather than an actual threshold) is not excluded by the slight 
inconsistency with the HI observations of the outer disc. 

4.3. The inverse chemical evolution problem 

The model of galaxy evolution described in Section|2]has many 
essentially free parameters. Nonetheless, models of this type 
have proven successful in several ways, e.g., by being able 
to simultaneously reproduce observational constraints such as 



the G-dwarf distribution, basic abundance patterns and the age- 
metallicity relation for stars in the solar neighbourhood, al- 
though the latter must be rega rded as a rather uncertain con- 
straint (see lFeltzing et al. 200 lb . However, due to the many pa- 
rameters, it is not obvious that one may take abundances derived 
from observations and, via CEMs, derive what stellar sources 
(types of stars) contribute the carbon and how the production of 
carbon in these stars may vary with metallicity. In principle, one 
has a really non-trivial inverse problem to solve, i.e., finding the 
stellar yields from observed abundance patterns. Basic require- 
ments for such an attempt are that a radically different star for- 
mation history should not alter the evolution of abundance ratios 
too much and that the IMF is universal and non-evolving or that 
it is possible to specify how the IMF is evolving. 

The modified sets of yields (C and D) were designed to im- 
prove the agreement between models and observations at dif- 
ferent metallicities. However, for stellar abundances in the solar 
neighbourhood, low metallicity often means that the star belongs 
to the halo/thick disc population. The halo/thick disc population 
has, of course, a different formation history and therefore the 
evolution of elemental ratios may be slightly different. Several, 
very diff'erent, infall-time-scales and star formation prescriptions 
were tested, but such variations of the CEM did not change 
the predicted solar neighbourhood abundance pattern much. The 
two different prescriptions for star formation used in this paper 
(type 1 and 2) provide very similar abundance trends for the so- 
lar neighbourhood with only minimal changes to other model 
parameters. In fact, the IMF and the stellar yields, or the IMF- 
weighted integrated yields, turn out to be the most critical in- 
gredients. Changes in the integrated yields will dominate over 
reasonable variations of all other parameters. It is therefore pos- 
sible to use CEMs to analyse how stellar yields may vary with 
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Fig. 9. Carbon abundance rel- 
ative to oxygen. Note that 
Galactic and extragalactic Hll- 
regions appear to follow a trend 
similar to that of the solar 
neighbourhood. The different 
lines have the same meaning as 
in Fig. 4-8, i.e., they correspond 
to models Al-El. The yellow- 
shaded (bright) area indicates 
roughly where C stars domi- 
nate the carbon enrichment of 
the ISM according to models 
Al, CI and El. Thin squares 
show data from Fabbian et al. 
(2009a) assuming no hydro- 
gen collisions in their abun- 
dance derivations, while thick 
squares show the same data 
set, but accounting for hydro- 
gen collisions. The quantitative 
model results are uncertain dur- 
ing the first few time steps, 
which is the reason the plot is 
cut at log(0/H) = 5.7 even 
though the data by Fabbian et 
al. (2009a) includes a couple 
of stars at even lower metallic- 
ity. The uncertainty is a conse- 
quence of the numerical meth- 
ods used and the fact that simu- 
lated abundance ratios may, in 
that early phase, change quite 
significantly due to the step 
sizes (linearly interpolated) in 
tabulated input data. 



Stellar mass and metallicity to be consistent with observations, 
provided that the IMF (and its possible evolution) is known. 

It seems that the carbon abundances in stars at low metallic- 
ity in the solar neighbourhood cannot be explained by a standard 
Galactic evolution scenario and commonly used stellar yields 
(cf. model Al and Bl). Consequently, modifying the yields to 
obtain better agreement with observations, as a solution to the 
inverse chemical evolution problem, may therefore be a justified 
measure. However, as model El shows, abandoning the idea of 
a universal and non-evolving IMF, may just as well be the so- 
lution to the carbon-enrichment problem for the early stages of 
Galactic evolution. Observed abundances (except those of ex- 
tremely metal-poor stars, perhaps) reflect the enrichment from a 
stellar population, and therefore one must be careful suggesting 
that abundances derived from observations pose a challenge to 
nucleosynthesis models and stellar evolution theory. We are left 
with two main options: (a) the carbon yields of essentially metal- 
free HM-stars must be much higher than predicted by WW95, or 
(b) the IMF evolves with time and/or is not completely universal. 
A necessary prerequisite for actually solving the carbon enrich- 
ment problem is thus to establish which hypothesis holds. In case 
both (a) and (b) is true, this would require a much less top-heavy 
IMF at early times, which may not be consistent with observa- 
tional constraints of some extremely metal-poor halo stars. 

The carbon enrichment at later stages (thin disc evolution) 
can also be explained in two ways, either: (a) the carbon con- 



tribution from HM-stars is metallicity-dependent in such a way 
that it increases significantly with metallicity during the thin disc 
evolution, or (b) the rising trend is due to a significant contri- 
bution from the long-lived stars (with masses < 2Mq) that be- 
comes important only after a few Gyr Scenario (a) is the one 
advocated by Gustafsson et al. (1999), w hile (b) represents the 
picture favoured by s everal authors (e.g.. IChiappini et al. 20031: 
lAkerman et al. 2004t iGavilan et al. 20051) . in which up to 80% 
of the carbon in the present-day ISM may be due to LIM-stars 
(see Fig.|5]l. It is virtually impossible to distinguish between (a) 
and (b) just by considering abundance data, but the fact that car- 
bon and iron seem to be released to the ISM on similar time- 
scales, suggests that these two elements originate from stars that 
evolve on similar time-scales. If LIM-stars produce most of the 
carbon and supernovae Type la produce the bulk of iron in the 
Galaxy, the rathe r flat trend of C/Fe with Fe/H would be ex- 
plained naturally dBensbv & Feltzing"2 006). If, instead, scenario 
(a) is correct, it would require that the metallicity dependence 
of the carbon yields of HM-stars is fine-tuned with the SNIa- 
produced Fe in such a way that both the flat C/Fe vs. Fe/H trend 
and the rising C/0 vs. O/H trend are reproduced. Thus, while 
(a) is a plausible hypothesis, scenario (b) is simpler in that it 
does not require that the yields have any specific metallicity de- 
pendence; it only requires that LIM-stars produce most of the 
carbon, so that the carbon enrichment of the ISM is delayed. 
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An interesting fact, which deserves further study (beyond 
the scope of this paper), is evident from data plotted in Fig. |9] 
Extragalactic Hll-regions show C/O-ratios which line up along 
the same trend with O/H as the stars in the solar neighbourhood: 
metal-poor (low-mass) galaxies show low carbon abundances, 
while more metal-rich galaxies show considerably higher car- 
bon abundances, relative to oxygen. This may, in principle, 
be consistent with both scenario (a) and (b) just discussed 
for the thin disc. Carbon production in the Galactic disc ap- 
pears to be coupled with the age/metallicity of the stellar pop- 
ulation. Similarly, it seems that the carbon abundance in ex- 
tragalactic Hll-regions is coupled with the age of the domi- 
nant stellar population or the overall metallicity. However, if 
(b) is the correct scenario, metal-poor galaxies must be domi- 
nated by young stellar populations, so that the majority of LIM- 
stars have not yet evolved off the main sequence, thus locking 
carbon, without significantly enriching the ISM. Nevertheless, 
the light distribution of HII/BCG galaxies reveals the presence 
of a significant low surface-brightness p opulation of old stars 
dTeUes et al. 1997l:lKunth & Ostlin 2000l) . although it is diflicult 
to determine the mass fraction of such an underlying population. 
If scenario (b) can be confirmed, metal-poor galaxies are very 
likely to have formed most of their stars quite recently, since a 
significant old population would raise the C/O-ratio due to en- 
richment from LIM-stars. 

It is often suggested (see iKunth & OstUn 200d[ and refer- 
ences therein) that low-mass, metal-poor galaxies are deficient in 
metals due to significant galactic winds and may therefore have 
undergone several star-forming episodes in the past. Galactic 
winds may be responsible for maintaining low over-all metallic- 
ities, but cannot easily explain low C/O-ratios unless the winds 
are extremely selective (oxygen remains while carbon is ex- 
pelled from the galaxy) in scenario (b). The fundamental ques- 
tion is then whether LIM-stars are (or not) the main contributors 
to carbon nucleosynthesis in metal-poor galaxies, and if galac- 
tic winds may maintain low over-all metallicities even after a 
sequence of star-forming episodes. 

4.4. Uncertainties in iron production 

Several uncertainties have not been addressed in detail here. 
One, particularly important, is the production of iron, which can- 
not be well-constrained for two reasons: the location of the mass 
cut (that divides the part of the star that collapses in the rem- 
nant from that which is expelled) in supernova nucleosynthesis 
models is a free (although restricted) parameter, and the origin 
of supernovae Type la (SNIa) is not fully understood. 

The iron yields by WW95 show a dependence on metallic- 
ity (possibly an effect of how the mass cut was chosen), which 
may not be realistic. AtZ-Q the IMF-weighted WW95 yields 
show an O/Fe-ratio that is much lower than the ratios derived 
from observations in very metal-poor stars ([0/Fe]~ 0.8 when 
[Fe/H]< -3, see Fig. |6]l, which is the reason why the iron 
yields at Z = have been lowered to 1/5 of their original val- 
ues in the models presented in this paper. However, non-local 
thermal equilibrium (non-LTE) effects on oxygen lines may 
explain the O/Fe discrepancy. Including corrections for non- 
LTE effects, smaller [O/Fe] values are derived from observa- 
tions, so the [O/Fe] may actually agre e quite well with WW95 
(iFabbian et al. 2006t iFabbian et al."20 09b). 

The fact that the iron yields of HM-stars are rather uncer- 
tain, and models of the evolution of SNIa rates cannot be very 
well constrained from observations, makes the common practice 
of using iron as the reference element (representing the over-all 



metallicity) somewhat questionable when comparing data and 
models. Hence, it makes more sense to calibrate CEMs primar- 
ily against abundances relative to, e.g., oxygen or other common 
elements for which the enrichment essentially follows the star 
formation rate. Therefore, focus is on reproducing C/H, O/H and 
C/O, rather than C/Fe and O/Fe, in the present study. 

4.5. Uncertainties in carbon production: effects of carbon 
star mass loss 

There are many assumptions and physical prescriptions that may 
affect the nucleosynthesis in stellar evolution models. For LIM- 
stars the duration of the AGB phase and the number of ther- 
mal pulses is almost uniquely determined by the mass-loss rate. 
The evolution of the internal structure depends on the mass-loss 
rate as well (Blocker 1995), which in turn affects the fundamen- 
tal stellar parameters. Since the mass-loss rate depends on these 
stellar parameters, there will be a feedback, which means that 
the mass-loss prescription put into a stellar evolution model is 
critical, and it is indeed well-known that changing the mass- 
loss prescription can have profound effects on the yields of AGB 
stars (see, e.g., van den Hoek & Groenewegen 1997). More pre- 
cisely, the sum of carbon and nitrogen synthesised in LIM stars 
is largely controlled by the mass-loss rate on the AGB. 

In the work of Gavilan et al. (2005) it was assumed that stars 
of low metallicity have smaller total radii (yielding higher sur- 
face gravity) and thus less effective mass-loss, which increases 
the lifetime of the AGB phase and allows these stars to experi- 
ence several more dredge-up events where carbon is mixed into 
the envelope. At late stages, as metallicity increases, the car- 
bon production declines, and secondary nitrogen pro duction be- 
comes increasingly significant (see also iBuell 19971 for further 
details). But is it correct to assume that low metallicity means 
less effective mass-loss relative to that of solar metallicity stars? 

Mattsson et al. (2008) have shown that, according to theoret- 
ical models, the overall metallicity is not affecting dust-driven 
mass-loss from C-stars much at all. It is mainly the abundance 
of carbon, available for dust formation, that has a significant ef- 
fect. The dust-driven mass-loss takes place mainly during the 
late stages of evolution, where thermal pulses and dredge-up 
events dominate the evolution. There is no particular reason to 
assume that low metallicity C-stars will experience more thermal 
pulses and dredge-up events before the termination of the AGB. 
Furthermore, Mattsson et al. (2010a) present a state-of-the-art 
model of the late stages of evolution of a 2M0-star using a new 
detailed mass-loss prescription based on the results by Mattsson 
et al. (2010b) and a new set of low-temperature opacity coeffi- 
cients, in which effects of the abundances of carbon and oxygen 
were taken into account. This resulted in a significantly lower 
effective temperature and the development of a very pronounced 
so-called superwind, shortly after the star had become carbon 
rich, which terminated the C-star evolution after only few ther- 
mal pulses. The amount of carbon that can be dredged up and 
expelled by the stellar wind is therefore limited. No final word 
yet, but these recent findings seem to restrict the carbon produc- 
tion in LIM-stars by means of a self-regulating process, where 
the mass-loss rate increases with every dredge-up event associ- 
ated with each thermal pulse. 

5. Conclusions and final remarks 

CEMs for the Milky Way have been presented. Four different 
nucleosynthesis prescriptions were used, where two contain ad 
hoc modifications to meet the abundance trends of C/O vs. O/H 
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and C/Fe vs. Fe/H derived from observations. According to these 
CEMs one may conclude the following: 

- Carbon is being released to the ISM on time-scales compara- 
ble to that of iron, which is produced mainly by supernovae 
Type la. 

- An evolving IMF, being top-heavy (favouring the formation 
of HM-stars) during the early stages of Galactic evolution, 
can explain the C/O vs. O/H trend seen in the abundances 
derived by Fabbian et al. (2009a), without violating any other 
observational constraints. 

- The C/H-gradient (or the metallicity gradient in general) de- 
rived from observations in the Milky Way disc suggest there 
is a problem with simple CEMs and the hydrogen distribu- 
tion. Introducing a star formation threshold at Ec - 1MQ 
pc"^ can explain the shape of the hydrogen distribution, but 
such a model predicts a steep fall-off in all abundance gra- 
dients beyond a certain galactocentric distance and cannot 
explain the flattening of the C/H and Fe/H gradients of the 
outer disc traced by Cepheid observations. It is possible that 
radial gas flows in the outer disc in combination with a 
"low-efficiency mode" of star formation is responsible for 
the abundance gradients and the hydrogen distribution. 

- HM-star yields of carbon and oxygen probably have a metal- 
licity dependence, as argued by Gustafsson et al. (1999), but 
the observed increase with metallicity of the C/O-ratio seen 
in disc stars could just as well be due to a delayed release 
of carbon from LlM-stars. If the best-fit models are correct, 
the major source of carbon in the present-day Galaxy is the 
LlM-stars, providing as much as 80% of the carbon to the 
ISM. Although HM-stars as major carbon producers cannot 
be excluded, the LlM-star scenario provides a simpler ex- 
planation, which does not require that carbon yields depend 
strongly on metallicity. 

The origin of carbon remains an open question. The scenario 
advocated here, is that LlM-stars have relatively high carbon 
yields and that the winds from these stars dominate the carbon 
enrichment of the ISM. However, recent results by Mattsson et 
al. (2010a) suggest that since the mass-loss of a carbon star has 
a steep dependence on the amount of carbon that is dredged up 
to the surface of the star, the mass-loss rate increases with with 
every dredge-up event, which eventually leads to the termination 
of the AGB phase. Hence, the number of thermal pulses will be 
very limited due to this self-regulating mechanism and the car- 
bon yield can hardly become as high as in models Dl and D2. 

That LlM-stars may contribute most of the carbon in the 
solar neighbourhood has p reviously b een concluded by several 
other authors (e.g., Chiap pini et al. 20 03: Akerman et al. 2004; 
iGavilan et al. 20051) . The C/Fe vs. Fe/H trend seen in unevolved, 
metal-poor stars in the solar neighbourhood suggests that the 
contribution from HM-stars at low metallicity must be limited, 
unless these stars also have an enhanced iron production, and 
the contribution has to increase dramatically with metallicity to 
explain the sharp up-turn in C/O vs. O/H, seen in Fig.|9] If HM- 
stars were the main contributors, then the carbon yields of these 
stars (around and above solar metallicity) must be even larger 
than the oxygen yields. Maeder (1992) and Portinari et al. (1998) 
have computed models where such results were obtained, and 
more recent work seems to give similar results (see, e.g., Meynet 
& Maeder 2002; 2003; 2005). But models of HM-star evolution 
are also riddled by several uncertainties. 

Finally, returning to the connection between the cosmic car- 
bon abundance and carbon-based life, one may notice that in the 



outer parts of the Milky Way disc, where the peak of star for- 
mation was reached quite recently (or has not yet been reached), 
there has not been enough time for the bulk of LlM-stars (espe- 
cially those with main-sequence masses below ~ 2Mq) to evolve 
into mass-losing C-stars, or the metallicity is not large enough if 
the HM-star carbon production scenario is instead the correct 
one. Are the outer parts of the Galaxy therefore excluded as a 
possible environment for complex life to emerge, since there has 
been no significant carbon enrichment until quite recently? 
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